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Abstract 

The branching ratios of the measured decay Kl — > 7r + 7r~e + e - and of the still 
unmeasured decay K + — ► 7r + 7r°e + e~ arc calculated to next-to-leading order 
in Chiral Perturbation Theory (CHPT). Recent experimental results are used 
to determine two possible values of the combination (N[ 6 — Nn) of weak low- 
energy couplings (LECs) from the 0(p 4 ) chiral Lagrangian. The obtained values 
are compared to the predictions of theoretical approaches to weak counterterm 
couplings to distinguish between the two values. Using the favoured value of 
the combination (N[ 6 — Nn) and taking into account additional assumptions 
suggested by the considered models, one obtains the branching ratio of the se- 
cond decay as a function of the unknown combination (N[ A + 2N[ 5 ) of weak 
low-energy couplings. Finally, using values of the individual LECs derived from 
a particular model, one predicts the branching ratio of the K + decay. 
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1 Introduction 



During the last years, there has been a lot of theoretical and experimental interest 
in the decay of the Kl into a pair of charged pions and a pair of leptons. This 
interest focused on the decay width itself ||, ||, [|, || |(| |], || and on the possibility 
of constructing CP-violating observables [2|, ||, f|, [j], ||, |To|, [ll], |l2| as well as on 



other related topics 14}| . 



From the experimental analysis of the corresponding radiative decay, it was found 
that the decay amplitude consists of a bremsstrahlung component and a direct emis- 
sion part. The contribution due to bremsstrahlung is given via Low's theorem by 
the amplitude of the decay Kl — > ir + ir~ . This amplitude is mainly due to the K® 
admixture, which allows for this decay (indirect CP violation). As a consequence, 
the final state of the radiative decay can be found to be in CP-even as well as CP-odd 
configurations. Hence, in principle, there is interference between the CP-conserving 
parts of the direct emission amplitude and the CP-violating bremsstrahlung ampli- 
tude. But as long as the polarization of the on-shell photon is not measured, this 
interference is not accessible. This is the reason why one looks directly to the decay 
with a lepton pair, since the angle between the two planes spanned by the pions and 
leptons can be used to construct a CP-violating observable jjj, |2[ |!| [| [R]] . 

In this paper, I do not focus on the CP-violating aspects of this decay. I calculate 
the decay amplitude in CHPT up to C(p 4 ) and use the most recent available data 
from experiments [L5], [L6| (which were mostly dedicated to the study of possible 
CP-violating effects) to derive a value for the unknown combination (iV{g(/u) — N17) 
of low-energy couplings (LECs) from the weak 0(p 4 ) chiral Lagrangian. Until now, 
theoretical predictions can only be compared to the branching ratio over the entire 
phase space, which makes it impossible to extract a precise value for this combina- 
tion. Furthermore, it is not possible to determine this value unambiguously from 
experiment, therefore one has to turn to LEC models and their predictions for low- 
energy couplings to find the favoured value. 

Once the value of this particular combination is fixed, I use it as input together 
with additional assumptions about the weak LEC Nu for the second non-leptonic 
decay discussed in this paper: K+ — > 7r + 7r e + e . If we use new data from the cor- 
responding radiative decay K + — > 7r + 7r°7 [17|, we can give the magnetic amplitude 
of K + — > 7T + 7r°e + e - without any unknown parameter at 0(p 4 ) and it is possible 
to predict the branching ratio BR(K + — > 7r + 7r°e + e~) as a function of the unknown 
combination (N{ 4 (fi) + 2N{ 5 (fi)) of weak LECs. 
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2 Effective chiral Lagrangians 



Chiral Perturbation Theory [[T^, |l9| is the ideally suited framework to discuss these 
processes. It is the low-energy realization of the Standard Model respecting the 
approximate chiral symmetry of the light quark sector. In fact, the demand of in- 
variance under chiral rotations (in our case, these are SU(3) rotations) allows one to 
write down the most general effective Lagrangian of strong interactions amongst the 
light pseudoscalar meson octet. The approximate chiral symmetry SU(3) l x SU (3)# 
seems to be realized a la Nambu-Goldstone, which means that it is spontaneously 
broken to the well-known SU(3)y. The breakdown of the symmetry gives rise to 
eight almost massless would-be Goldstone bosons because there are eight broken 
axial generators. According to Goldstone's theorem, the quantum numbers of these 
particles are fixed by the quantum numbers of the broken generators, thus one iden- 
tifies the light pseudoscalars with these particles. 

In the scheme of Gasser and Leutwyler |l9|], the most general 0(p 2 ) Lagrangian 
including strong, electromagnetic and semileptonic weak interactions reads as fol- 
lows: 

4 

where D^U is the covariant derivative with respect to external, non-propagating 
fields. If we specialize to the case of external photons, 



C 2 = —(D tl UD^ + X U^ + X^U), (1) 



Q = i-diag(2,-l,-l), (2) 

where Q is the quark charge matrix for the flavours up, down, strange. U is a 3 x 3 
unitary matrix which has to be expanded to the relevant order in <I>: 

U($) = e <V2*/^ ( 3 ) 
where the mesons are collected in the matrix 



V2 



(4) 



F equals to lowest order the pion decay constant, F n = 92.4 MeV. In general, \ 
contains external scalar and pseudoscalar matrix-valued fields, but here it is propor- 
tional to the quark mass matrix. In this way, explicit chiral symmetry breaking can 
be incorporated in the effective Lagrangians in a very elegant way: 

X = 2B ■ diag(m u , m d ,m s ). (5) 
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-Bo is related to the order parameter of the spontaneous breakdown of the chiral 
symmetry, the quark condensate. It will not appear explicitly because it can be 
absorbed in the squared meson masses. 

For the calculation of non-leptonic kaon decays, we also need an effective La- 
grangian describing the weak interactions of the mesons. This effective weak La- 
grangian cannot be invariant under chiral rotations, hence chiral invariance cannot 
be the guideline. Starting from an effective strangeness-changing AS" = 1 four-quark 
Hamiltonian, one writes down a hadronically realized Lagrangian that transforms 
in the same way under SU(3)l x SU(3)r as this Hamiltonian. At lowest order, the 
needed weak Lagrangian is found to be: 

£* s=1 = G 8 {XL^)+G 27 [L^ 1 + h^ 1 L^] + h.c., (6) 

where A = (A6 — i\j)/2 projects out the correct octet quantum numbers and 
= iF 2 U^ D^U is the hadronic left-chiral current in analogy to the left-chiral 
quark current at the level of the effective Hamiltonian. The two couplings G$ and 
G27 have to be obtained from experiment and the determination of these couplings 
involves some subtleties. In principle, the couplings are obtained from K — > irir 
decays. Comparison of experiments with the leading order 0(p 2 ) calculations yields 
the 'canonical' values \G$\ ~ 9.1 • 10~ 6 GeV~ 2 and G 2 7/G$ ~ 1/18, where this 
approximate ratio of the two couplings introduces uncertainties when the 27-plet 
coupling enters the game. However, due to the smallness of the 27-plet coupling, 
one can usually neglect this part of the Lagrangian unless the octet contribution 
vanishes. Then also the 27-plet contribution may become important (see Sect. 3.2). 

For completeness, one should remark that in Ref. [pOf | the relevant K — > -kit 
decays were analyzed up to 0(p 4 ) and it was found that these additional corrections 
contribute to Gs with about 30%, whereas the G27 coupling is only modified by a 
few percent. Thus, if one takes into account these order p 4 corrections, the value of 
the coupling | Gs | appearing in (|) should better be ~ 6.4- 10" 6 GeV" 2 . Throughout 
this work, however, I am using the canonical standard values. 

The chiral Lagrangians (|l]) and (|6|) allow us to calculate tree-level amplitudes 
of chiral order p 2 and one-loop diagrams of chiral order p 4 which usually introduce 
divergences. In order to get rid of these divergences and to take into account further 
finite local corrections appearing at 0(p 4 ), e.g. through new interactions arising 
from the chiral anomaly, one also has to consider the most general 0(p 4 ) interaction 
Lagrangians. 

The most general strong Lagrangian of order p 4 , invariant under C, P and chiral 
transformations, was again given by Gasser and Leutwyler There is only one 
term in this Lagrangian that contributes to the final results in this work: 

£4 = -iWFgDJJDyrt + F^D^D U U). (7) 



3 



Since we are interested in external photons, the F^ U R tensors are proportional to the 
ordinary electromagnetic field strength tensor: 

Fjf = -eQF^ u = F£, F» v = d»A v - d u A^. (8) 

Every new term in the strong Lagrangian of order p 4 is furnished with an a priori 
unknown low-energy coupling (LEC) [|19| . Since all divergences appear as local 
polynomials, one can absorb the divergences of the loop amplitude in the LECs. 
The general structure of a LEC reads as 



d-4 r 



Li = LKjm) + T 4 A( M ), A(^) = — — - -(ln(4vr) + 1 - lE ) , (9) 



16vr 2 



1 



d-4 2 



where "Je = 0.5772157 is the Euler-Mascheroni constant. This is also true for weak 
LECs Ni . The coefficients T j arise from the one-loop generating functional. Because 
of the regularization procedure, the measurable couplings L\{^) (and N[(/j,)) become 
scale dependent. In the end, this scale dependence must be compensated by the scale 
dependent parts of loop diagrams. One should also note that the chiral subtraction 
prescription differs from the usual modified MS prescription. 

The new octet weak interactions are organized like this |21] , f22| : 

Cf S=1 = G$F 2 ]T NiWi + h.c. (10) 

i 

For the non-leptonic kaon decays under consideration, only the operators W14, W15, 
WiQ,Wn and W28, ^29, W30) W31 contribute; they are listed explicitly: 

Wi 4 = i{\{F£ , + U*Fg'U,D ll rt D V U}), 

W 15 = i(XD^(UF^U^ +F R W )D V U), 

W w = i(\{F£ u — U^F^UjDfJJ^DjjU}), 

W 17 = i(XD^(UFfU^ -Fg)D v V). (11) 

The magnetic terms (proportional to e pupa ) are given by 

W 28 = ie pupa (\D^U)(U ] D U UDPU ] D a U), 

W 29 = 2(X[U j 'F R M U,D^D u U]}, 

W 30 = {XrtDpU)^ + rfFgU)D v rtU), 

W31 = (\U^D p U)((F^ -U^F^ V U)D V U^U), (12) 

with F% R the dual tensor of (|), F^ R = F paLiR . 

Finally, we introduce a Lagrangian that embodies contributions from reducible 
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diagrams with a WZW-vertex and an 0(p 2 ) AS = 1 vertex. It only contributes to 
the K + decay and is given by H^, [24], |25| 

£t S=l = ^F^dy K +15„n-, (13) 

where F^ v is the dual of the ordinary electromagnetic field strength tensor (||), 
pvv _ e ^v P ap an( j ^ e covar iant derivative is the usual QED derivative. 



3 Amplitudes 

For both decays, the general form of the invariant amplitude due to covariance is 

A = ^V lt u(k-)-fv(k+), (14) 

where q = + k + is the momentum of the virtual photon, /c_ and k + are the 
momenta of the electron and positron, respectively. iV^ is the generic weak Kmr^*) 
vertex, calculated in CHPT. It is decomposed in an electric and a magnetic part 

Vfi = FiPip + ^2P2 M + Me^paplp^q ', (15) 

where pi andj>2 are the outgoing momenta of the ir + and 7r _ (7r°) and T\, J- 2, M are 
form factors containing the dynamics of the two processes. A separate term propor- 
tional to the photon momentum vanishes because of the Dirac equation. The form 
factors are either constants or scalar functions of various products of the involved 
momenta. 



3.1 l\ 1 — > 7T+7T 7* Amplitudes 

This decay had already been considered in the framework of CHPT in ||. The 
authors of Ref. || used a different basis of counterterms (this change of basis 
is only valid as long as one is only interested in photons in @) and a different 
approximation of the magnetic part of the amplitude not taking into account any 
energy dependence. The present calculation considers this energy dependent part 



11 ], too, and additionally serves as a check on the results in Ref. H. 

For this decay, I assume strong isospin conservation, i.e. the up and the down 
quark have equal masses. Hence, the Gell-Mann-Okubo mass relation holds and we 
only have to deal with two independent masses: 3m^ g = Am^ — m 2 . It will be used 
to simplify parts of the one-loop amplitude given in the appendix. 

In this paper, I use the following definitions: Kl = K® + eKf , where CP\K®) = 
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K°(p) 




Figure 1: Tree-level diagrams for Kl 
entirely due to Kf admixture. 



7T+7]- j*. At tree level, the Kl transition is 



+ \K®) and CPlK®} = —{K^). K® and K\ are related to the strangeness eigenstates 
K° and K° through the following expressions: 



K 



Kn 



(16) 



The tree-level amplitude is entirely due to the admixture, since we do not 
consider direct sources of CP violation. In any case, the tree-level contribution is 
rather suppressed, especially when compared to the K + decay (Sect. 3.2). From 
Figure [l] one obtains the following tree-level form factors: 



Lt 



-ie 



4eG 8 F 
2qpi + q 2 



[m 2 K 



Lt 



. JeGsF 2 
ie- -{m K 



(17) 



where e ~ 2.27- 10 _3 e* 44 ° is the parameter of indirect CP violation. In the remainder 
of the paper, we do not take into account 0(p 4 ) corrections proportional to e to the 
electric form factors. As mentioned in Sect. 2, a value of |Gg| ~ 9.1 • 1CP 6 GeV -2 
already amounts to some 0(p 4 ) contributions. 

The magnetic form factor can only arise through the four weak counterterms W28, 
. . . ,Ws\ and it is in fact a result of the chiral anomaly. It is necessarily finite and 
does not have any energy dependence at this order. CHPT generates the following 
direct emission magnetic form factor: 



M 1 



= ^2^( a 2 + 2a 4), 



(18) 



where I have used the 'magnetic' notation of Refs. [24, 25]. These magnetic LECs 
are also still unknown. Experiments exhibit a large sensitivity of the magnetic 
amplitude to the energy of the emitted photon, therefore I will use the experimental 
results of [11] (rather than the old results of the experiment by Ramberg et al., [p^ ]) 
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to estimate the magnetic contribution. The authors of Ref. [11] use the papers by 
Sehgal et al. (I], ||] as the theoretical background to model their Monte Carlo, but 
additionally introduce an energy dependence in the magnetic amplitude through a 
form factor that involves a kind of a p propagator: 



,L „ UM 5Afl- 



m 



K 



ai/a 2 



M = e|/ s |^W, W = 1 + — 9 . (19) 



(m 2 p - ra\) + 2m K E* 



Ansatz (19) cannot be compared directly to the magnetic form factor in [|], ^, ||, fl]. [f| . 
Consequently according to 11], one should identify the average of <7miVV over the 
allowed range of E*, the energy of the virtual photon, with the original mag- 
netic coupling used in Refs. 0, g, |, g, g|. |/,| ~ 3.9 • 10~ 4 MeV is the abso- 
lute value of the decay amplitude of Kg — * ir + ir~, and the experiment [11] gave 
for the magnetic coupling \(Jmi\ = 1.35lQ^(stat.) ± 0.04(syst.) and for a\ja 2 = 
-0.720 ± 0.028(stat.) ± 0.009(syst.) GeV 2 . (a 2 in the fraction above is not the 
same as the LEC a 2 in (|l8|).) These numbers were obtained from the entire KTeV 



1997 data set of more than 1811 events above background [11]. In fact, it is also 



this data set and this parametrization that were used to extract the most recent 



value of the branching ratio of Kl — > 7r + 7r~e + e~ ||15| , |16f l- Additionally, the frac- 
tion a\/a 2 was found from the corresponding radiative decay Kl — ► 7r + 7r _ 7 to be 
— 0.729 ±0.026(stat.) ±0.015(syst.) GeV 2 [12|, which is clearly in perfect agreement. 



The errors of these quantities are the sources of by far the most important contribu- 
tions to the uncertainties in the extraction of the LEC combination (A r [ 6 (^) — Nn). 

The electric form factors at 0{p A ) show the pleasant feature that one can obtain 
the form factor J 7 ^ fr° m the expression for by simply exchanging the pion mo- 
menta pi and p 2 . At this order, there is no change of sign as at the tree level (|T^), 
since we concentrate on the CP-conserving part of the decay. 

We begin the discussion of next-to-leading-order electric form factors by con- 
sidering strong loops and strong counterterm contributions. The starting point for 
our analysis is the collection of diagrams in Figure |], where one replaces with 
K\- Wave function renormalization graphs will be neglected, since the tree-level 
amplitude of K® — ► 7r + 7r~7* vanishes. 

Removing the photon in the left diagram of Figure [l] and replacing the external 
kaon line with the appropriate loop diagrams from Figure |2|a, one obtains diagrams 
that are found to vanish. This feature is due to the structure of the weak i^2 7r+7r_ 
vertex. Appropriate replacement of the kaon line in Figure |l] with the loop diagrams 
drawn in Figure [| however, yields a finite propagator contribution to the electric 
form factors. The contribution must be finite, since there are no counterterms to 
compensate a divergence. This kind of diagrams was already considered in [Q, |2| 
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(a) (b) 

Figure 2: Strong insertions, a) Loops, where the photon is emitted at the vertex 
or by a charged meson in the loop, and the generic photon emitting counterterm 
proportional to Lj. M denotes K®, vr + or 7r~. N denotes the allowed particles 
in the loop: 7r°, 7r + , K + , r/g, K® or K\- O denotes any charged pseudoscalar. P 
denotes ir + or ir~ . b) Generic loop and generic counterterm vertex proportional to 
Li without a photon. 




Figure 3: Strong loop insertions for the external K% line that generate a propaga- 
tor. The contributions derived from these insertions, however, are finite. O denotes 
a charged pion or kaon. 



and the sum of the diagrams yields 

(d - 1)F [(pi +p 2 y - m z K \ 
{B(q 2 ,ml,ml)(Aml - q 2 ) + (4 - 2d)A{va\ ) 

-B(q 2 ,m 2 K ,m 2 K )(Am 2 K - q 2 ) - (4 - 2d)A(m 2 K )). (20) 

A second kind of strong 0(p 4 ) corrections is obtained from the two remaining 
bremsstrahlung graphs in Figure |l] in two ways: first, by either putting the loop 
diagram (without photon) of Figure |2|b or the counterterm insertion (without pho- 
ton) of Figure Qb instead of the internal pion lines. Secondly, by replacing the scalar 
QED vertices with the diagrams of Figure ^a. Focusing on the counterterms first, 
one finds that counterterms proportional to the LECs L4, L5 and Lg from the order 



p Lagrangian [19| are allowed to contribute. Calculating the sum of all the dia- 
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grams of the second kind, however, one discovers that only the contribution from 
the counterterm (0) proportional to Lg survives. The final correction of strong order 
p 4 diagrams of the second kind to the form factors is given by a very condensed and 
compact result: 



J~ 12 



—ieG% 



F 



{ - Vl, 



d-2 



+ 



{Ami 



2A(ml)+A(m 2 K ) 
q z )B(q z ,ml,ml) 



d-l 

,2\ or „2 2 



+ ^m 2 K -q 2 )B(q 2 ,m 2 K ,m 2 K )]} 



(21) 



The tadpole integral A(m 2 ) and B (p 2 , m 2 , m 2 ) , the scalar two-propagator integral, 
are defined in Appendix A. d is the spacetime dimension coming from dimensional 
regularization. 

Actually, the form factor contains divergences stemming from the functions 
A and B which are removed by the divergent part of the strong counterterm coupling 
Lg (coming from Figure |^a). In the finite amplitude with strong 0{p A ) insertions, 
the measurable part of Lg shows up: L r 9 {m p ) = (6.9 ± 0.7) • 10~ 3 Jl9| ]. Throughout 
this paper, I always choose fi = m p as renormalization scale. 

Weak counterterms only contribute through a diagram obtained from the direct 
emission diagram in Figure [l] by replacing K® with K® and putting in the coun- 



terterm vertex from (10) and ( |ll|) instead of the lowest-order vertex. All occuring 
divergences from weak loop diagrams must be removed by this local counterterm 
contribution. The weak counterterms produce the following contribution to the 
electric form factors: 



■r-Ll 

■r 13 



2ieG. 



8 2 



3F 



[N U -N 15 -3(N 16 -N 17 )}, 



(22) 



rn 



where all renormalized LECs N[ (compare with (||)) depend on a scale \i 
Since the coefficient Tyj is found to vanish, the LEC iVi7 is independent of the 
renormalization scale. The renormalized parts of the low-energy couplings enter 
into the amplitude of the decay, therefore it is important to know their finite values. 

The combination (A^{ 4 — N^ 5 ) also appears in the counterterm part of the form 
factor describing the decay K + — > ir + e + e~ within the expression [27] 



64tH 



w + 



[N[M-N r M + 3LM] + lln 



(23) 



Old experiments [^] fixed w + to be 0.89^^14, which corresponds to a value of 
{N[ 4 {m p ) — Nl 5 (m p )) ~ —0.02. A more refined theoretical analysis of this decay 
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also took into account 0(p ) corrections to the form factor |29[| . The polynomial 
part of this form factor is given by W+° l = GFm 2 K (a + +b + z), where z = q 2 /mn? K and 
q is the momentum of the intermediate photon that gives rise to the lepton pair. 
The new parameter a+ contains in principle also 0(j> 6 ) corrections and it is related 



with the usual w + through [29] 



r/.+ 



Gf 



w+ 



(24) 



A new experimental analysis of this decay [3C] measured the parameters of the 
K + -> vr+e+e" form factor and found a+ = -0.587 ± 0.010. With this new number 



we determine w + to be 1.086 and (iVj" 4 (» - N{ 5 ([i)) = -0.019 ± 0.002 at the scale 
m p . One finds that the new and the old value are almost the same. 

The contributions of weak loop graphs to the form factors are quite involved. 
To make it more transparent how the corrections from different kinds of weak loop 
graphs enter into the form factors, I present the possible kinds of diagrams in Figures 
[| and |5| and quote the results separately. 




Figure 4: Weak loop diagrams: the basic tadpole diagram (left) and the basic 
diagram of topology 1 (right). N and O denote the same particles as in Figure ||. 

All weak tadpole diagrams can be obtained from the basic diagram (left) in 
Figure ||| by appending a photon on all charged lines and on the weak vertex. It 
turns out that only intermediate charged particles produce non-vanishing diagrams. 
The tadpole part of the form factor looks very simple and reads as 

= = y 8 d ~ l {2A(m£)(2d - 4) + 2B(q 2 ,ml,ml)( q 2 - Ami) 

+A(m 2 K )(2d - 4) + B(q 2 , m 2 K ,m 2 K )(q 2 - Am 2 K )}. (25) 

Diagrams which can be constructed from the second diagram (right) in Figure || are 
referred to as diagrams of topology 1. Again, one has to append a photon on all 
charged lines as well as on the strong and weak vertex. This time, only pairs of 
charged pions or charged kaons may occur in the loop. The contribution of topology 
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1 to the form factors is found to be very compact, too, and it is given by: 



- d)A(ml) - B(q 2 ,ml, ml)(q 2 - Ami) 

±B(q 2 ,m 2 K ,m 2 K )(q 2 -4m 2 K )}. (26) 

The diagrams considered so far produce form factors that are symmetric in the pion 
momenta p\ and p2- Besides, apart from A functions only B(q 2 ,m 2 ,m 2 ) occurs and 
one can easily check that all these contributions vanish for an on-shell photon. 

The decay amplitude is completed with the contributions from diagrams be- 
longing to topologies 2 and 3. These diagrams are obtained from the basic graphs 
in Figure || through the same steps as before. The expressions that one obtains 
from these graphs are rather involved, thus I will not present the results explicitly 
in terms of the standard scalar loop functions A, B, C defined in Appendix A. It 



u -ieG 8 

>15 " 



^{2(2- 
+ (2 - d)A{m 2 K ) - ■ 




Figure 5: Weak loop diagrams: the basic diagram of topology 2 (left), the basic 
diagram of topology 3 (right). 

is the contributions of these diagrams that introduce the asymmetry in p\ and p2 
in the C(p 4 ) form factors. The possible pairs of particles in the loop are (jr ,K~), 
(r]8,K~), (Kf,Tr~) and (K%,it~) for topology 2. The particles for topology 3 are 
the corresponding charge conjugated ones. It turns out that the diagrams with the 
internal combination (K^ir') vanish. The contributions from the other possible 
combinations to the form factor T\ are given in J-^ in Appendix B, Eq. (f4^). 

By extracting only the explicit poles of the total loop contribution, one finds 
that all divergences are proportional to q 2 , which corresponds to the counterterm 
parts of expressions (|2l|) and fl22]). Furthermore, this shows that the loop amplitude 



of the corresponding radiative decay is finite [23, 25]. Finally, the complete form 
factor T\ is given by 

T{ = T{ 1 + T[\ + T{\ + F{£ + T[\ + T[l + T{1 . (27) 

As already stated, the corresponding form factor T2 is obtained from T\ through 
the substitution p\ <-> p2- 
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3.2 K — > 7r+7r°7* Amplitudes 

The general structure of the amplitude stays the same as in 

(U), but there is no 

symmetry relation between the electric form factors anymore. p\ is now the momen- 
tum of the 7r + and P2 the momentum of the 7r°, respectively. In the limit of isospin 
symmetry, the octet tree-level amplitude vanishes, hence we relax the approximation 
of equal masses of charged and neutral pions at the tree level and take the 27-plet 
coupling into account, too. As already anticipated in Sect. 2, throughout the fol- 
lowing analysis we will use again the canonical values of | Gg | and G27 /G% which are 
derived from the tree level. 

The tree-level form factors arise from the corresponding diagrams in Figure [l], 
where one replaces K® with K + and ir~ with tt° and puts the photon into the 
right places. In addition, the corresponding tree-level contributions with the weak 
coupling constant G27 will be regarded, too. Terms proportional to Gs are clearly 
suppressed because of approximate isospin symmetry, thus the actual value of Gg is 
not of too much importance for the tree level. The lowest-order amplitude reads as 

J? = «eGaF(^-^){^ r? + ? ^} 

+ ^^(5< + -7^ ++ 2^o), 



3 V K+ ^ ^ J (qZ + 2q Pl )(ql-2qpy 



2ieF 



2 ~~ n2 



q 2 — 2qp 



G 8 (ml + -ml,)- ^-f- (5m 2 K+ - 7ml + + 2m\ )}. (28) 



The magnetic form factor at lowest order (C(p 4 )) is derived from the Lagrangian 
fllOD with the counterterms in ([E|) and the WZW Lagrangian (|l5|). It is necessarily 
finite and one calculates 

• M+ = ^l (2 " 3a2+6a3) ' (29) 

where the 2 comes from the Lagrangian (|l^). Again, the values of these magnetic 
LECs are unknown, but comparison with the corresponding radiative K + decay 



24, Eq] shows that the magnetic form factor A4 + also appears there. This suggests 



to use results from the E787 experiment [17] to estimate the combination of LECs in 



(29). In this experiment on the corresponding radiative K + decay, a branching ratio 
from direct emission BR(K + -» vr + vr 7; DE, 55 MeV < T 7r + < 90 MeV) = [4.7 ± 
0.8(stat.) ± 0.3(syst.)] • 10 -6 is reported. Under the rather reasonable assumption 
that direct emission is entirely due to the magnetic amplitude, one can extract a 
value for the whole combination of LECs in (|29|). Of course, this does not take 
into account energy dependent corrections, but this is at the moment the best one 
can do. Moreover, the experimental data seem to indicate that neglect of energy 
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dependent higher-order terms does not do much harm to the magnetic amplitude. 



The authors also find no evidence for any electric direct emission in the decay [17|. 
The combination of magnetic LECs in (p9[) can be extracted from the radiative decay 



(q 2 = 0) by using A(K + 



7r + 7r°7, DE) 
3a 2 + 6a 3 | = 



= M<< 
\AA = 2.26 ±0.25. 



(30) 



Turning to next-to- leading-order corrections to the electric form factors, I start 
again with the discussion of contributions from strong loops and strong counterterm 
diagrams. From now on, strong isospin is conserved and 27-plet corrections are ne- 
glected. 

Similar to the analysis in Sect. 3.1, diagrams with strong loops and strong coun- 
terterm vertices are obtained from insertions of loops and vertices in propagators or 
external lines in K + — > 7r + 7r°7* tree- level diagrams derived from the graphs in Figure 
|l] by replacing Kf — > K + and 7r _ — > 7r°. Wave function renormalization diagrams are 
not considered because the tree-level octet amplitude for K + — * 7r + 7r°7* vanishes in 
the isospin limit. The necessary insertions are obtained from the diagrams in Figure 



and N denotes ir + , ir , K+, 



or 



where M denotes this time K + , tt + or n 
t/8, respectively. Here, P denotes K + or ir + . Introducing equivalent replacements 
as in Sect. 3.1 one obtains the strong corrections to the next-to- leading-order form 
factors given by: 

—ieG» r . o _ d — 2 r 



•77+' 

•> 1 
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+ 



F 
1 

d-l 



Aq 2 L 9 



d-l 



2A(mi) + A(m z K ] 



(Ami ~ <?)Btf 



,m n ,m nj 



+ \{Am 2 K 



,m 2 K ,m 2 K , 



)]}• 



T 

J r 



+1 



21 



0. 



(31) 

The vanishing of is related to C invariance of the strong Lagrangian. Although 
counterterms proportional to L4, L5 and Lg are allowed to contribute, only the Lg 
term (^) survives the summation of all contributions. The divergent part of the LEC 
Lg removes the divergences of the result in (31). The form factor J 7 ^ 1 vanishes in 
the limit of an on-shell photon. 

Substituting the relevant weak counterterm vertex from ( |l0| ) and (|ll]) for the 
lowest-order vertex in the direct emission diagram in Figure |l] and making the neces- 
sary particle replacements, one calculates this local contribution to the form factors: 



f+l _ 

•M2 — 


—ieGg 
3F 


T +l _ 
x 22 — 


—ieGg 


3F 



N- 



- 6qP2(N u - N 15 - N 16 - N 17 ) - Aq 2 (N 1A 

6q Pl (N u - N 1B - N w - N 17 ) - 2q 2 (N 14 + 2N 15 ) 
+6q 2 (N m -N 17 ) . 



15, 



(32) 
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One recovers the finite combination (-ZV14 — A15 — ^16 — A17) of the corresponding 
radiative decay and the structure that is governed by gauge invariance 25 1. 
Divergences arising from weak loop diagrams are removed by the combinations of 
LECs proportional to q 2 in ([32]). Using the determined value of (N[ e — Nn) and 
appealing to some models for weak low-energy couplings, the whole finite decay 
amplitude contains in the end only the unknown combination (N'( A + 2iV{ 5 ) from the 
form factor J 7 ^ m fl3^). I will come back to this later. At this point, it should be 
mentioned that a similar combination of the same weak LECs, namely (2-/VJ 4 + N[ 5 ), 
appears in the decay Kl — ► 7r°7r°7* considered in Ref. [31]. 

Weak tadpole diagrams can be constructed from the basic diagram in Figure |I| 
(with K\ — ► K + , ir~ — > 7r°) by following the same procedure as in Sect. 3.1 and one 
obtains 

HI = "3^ 8 / 1 {(2d-4)[AK)+A(m^)] 

+B(q 2 ,ml,ml)(q 2 - im 2 n ) + B(q 2 , m 2 K , m 2 K ){q 2 - 4m 2 K )}, 

+B(q 2 ,ml,ml)(q 2 - Ami) ~ B(q 2 ,m 2 K ,m 2 K )(q 2 - Am 2 K )}. (33) 

One finds that rjg loops do not contribute. In case of an on-shell photon, expressions 
( p3| ) vanish. 

Diagrams of the topology 1 are constructed from the right diagram in Figure 
U (with K\ — ► K + , ir~ — > 7r°) replacing the charged meson pairs in the loop with 
(vr + , 7T ), (i^ + , -K -0 ), or with (7r + , r/s) (in an appropriate momentum convention). The 
last combination of intermediate particles vanishes in the isospin limit. Appending 
a photon where it is possible and summing up the diagrams, one finds a compact 
result involving only A and B(q 2 ,m 2 ,m 2 ): 

ieGx r , , „,r 9 s 4 



HI = j^T ) {(d-2)[A( m 2 K ) + U( m 2 7r ) 



+ ^B(q 2 ,m 2 K ,m 2 K ){q 2 - 4m 2 K ) + ^B(q 2 , m 2 , ml ){q 2 - Am 2 )), 

Hi = j^^{(2-d){A( m 2 K ) + lA( m l) 

+ ^B(q 2 ,m 2 K ,m 2 K ){4m 2 K - q 2 ) + ^B(q 2 ,ml,ml)(4ml - q 2 )}. (34) 

Expressions (|34| ) vanish in the on-shell limit and their divergences are clearly pro- 
portional to q 2 . The last and by far most voluminous contributions to the electric 
0(p 4 ) form factors come from diagrams of the topologies 2 and 3 which can be 
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derived from the basic diagrams in Figure || as before. Possible virtual pairs are 
(tt°,K-), (t] 8 ,K-), (K°,ir-) and (K°,ir°), (K°,r] 8 ), (K+,ir + ), respectively. The 
obtained results, labelled as an d ^16 26> are listed in Appendix C, 

and (|5^), (|5l|). The complete form factor J 7 ^ is finally given by 



F+ = Tt l + F# + ^l + 2 + ^l + 3 + ^l + 4 + ^l + 5 + ^l + 6 • (35) 

J- 2 is obtained from the corresponding sum. 

4 Numerical analysis 
4.1 Decay width 

The decay width for the processes in question is given by the following standard 
formula: 

„ fir + m 2 e f d 3 p l d 3 p 2 d 3 k + d 3 k^ (4) 

~ spins 

(36) 

where pi is always the momentum of the positive pion and ^2 refers to the cor- 
responding other pion, 7r~ or ir°. As usual, pjj denote the sums of ingoing and 
outgoing momenta, respectively. In fact, pf = (m^, 0, 0, 0). The squared transition 
amplitude for both decays in question reads as follows: 

^2 



El" 4 !' = ^ T {-K + ^+^)fl^llV + l^2| 2 ^+PlP2(^+^2) 

+2\F 1 \ 2 k +Pl k_p 1 + 2\T 2 \ 2 k + p 2 k^p 2 + (TxFl + T{T 2 ) 



spins 



2 2 
Pl<7£>2 



1 e f 
(k + pik-p 2 + k-pik + p 2 ) \ H 5-^-^ (-mj + /c+fe_ 

+P2QP1 + Q 2 PlP2 ~ PipW - 2piP2qPiQP2 + 2qk-qk + 
(P1P2-P1P2) + 2k -Pik + pi{plq 2 - qp\) + 2k^p 1 k + p 2 (pjq 2 - qp\) 
+2(qp 1 qp 2 - q 2 pip 2 )(k + p 1 k_p 2 + k + p 2 k^px) + 2(p 1 p 2 qp 2 - p\qp\ 

(k + pik-q + k-pxk+q) + 2{pip 2 qp 1 - p\qp 2 )(k + p 2 k^q + k_p 2 k + q) 



+ 



m^q 4 



^ pa k_ uPluP2p k +a {(k +Pl - k- Pl )(J^M + FiM*) 



+{k +P2 - k_p 2 ){F 2 M + F 2 M*)\. (37) 
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The structure of (37) implies that there is no interference between electric and 
magnetic form factors in the decay widths of these decays. Additionally, one finds 
for the decay of the Kl that there is no interference between electric form factors of 
lowest and next-to-leading order, too. This feature is due to their different behaviour 
under exchange of pion momenta. Kl branching ratios thus consist of three distinct 
contributions. The more general case of interference between electric form factors of 
different orders is present in the K + decay. Phase space integrations are performed 
numerically with the Fortran event generator RAMBO [B2I|. 



4.2 Numerical analysis of Kl — > n + ir e + e 

In the following, branching ratios (BRs) for different cuts in q 2 , i.e. for different 
lower bounds on (£;+ + k-) 2 , and the BR over the entire phase space are listed, q 2 
may vary between Am 2 and (mx — Iva-^ 2 . Throughout this analysis, the central 
values of experimental numbers are used for the branching ratios with certain cuts 
in q 2 . The error of a branching ratio is only given if it is calculated over the whole 
phase space. It should be pointed out (compared to Ref. [||]) that the KTeV data, 
on which I will rely in the following analysis, are corrected for the entire phase space 

PL II, El- 



Table 1: Magnetic and tree-level contributions to the branching ratio of Kl 
7r + 7r~e + e~ for different cuts in q 2 and for the entire phase space. 

q 2 > (MeV 2 ) Magnetic BR [1(T 8 ] Tree-level BR [1(T 8 ] 



2 2 


18.20 


9.8 


10 2 


9.31 


2.95 


20 2 


5.61 


1.33 


30 2 


3.65 


0.71 


40 2 


2.42 


0.41 


60 2 


1.06 


0.16 


80 2 


0.44 


0.061 


100 2 


0.16 


0.024 


120 2 


0.053 


0.009 


180 2 


0.00025 


0.0001 


entire p.s. 


21.2 ±9.0 


12.8 ± 1.0 



Using ansatz ( |l^) and the experimental numbers of jTl]] as input for the magnetic 
contribution to the branching ratio, I find the results collected in Table M. Neglecting 
the energy dependent part in (|l^) one reproduces the results in Ref. || . Obviously, 
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consideration of the energy dependent magnetic form factor in ( |19| ) increases the 
results compared to a constant magnetic form factor, particularly for low cuts Q. 
Unfortunately, the errors of the parameters entering into the magnetic contribution 
to the BR are rather large [11|, thus the magnetic branching ratio over the whole 
phase space has a considerable uncertainty. 

The tree-level form factors of ( |l7|) give rise to the results collected in column 
three of Table |]. Comparison with the results in || shows that the obtained numbers 
are rather different, but this is due to different values of F and Gs- Here, we use 
F = 92.4 MeV and the canonical \Gs\- The error of the BR over the entire phase 
space in the last line comes from numerics and reflects the l/q 4 behaviour of the 
squared amplitude. Table [l] also shows very clearly the importance of the q 2 range 
between 4m 2 : and 4 MeV 2 that was not considered in The importance of this 
small q 2 range is understood from the plot of the differential decay widths of the 
individual contributions to the decay in Figure JBL 

The O (p 4 ) electric form factors depend via (|^) on (iV{ 4 (/i) — N[ 5 (fi) — 3(N[ e (//) — 
Nn)) =: X(fi), hence electric next-todeading-order contributions to the branching 
ratios are given as functions of X. The derived branching ratios are listed in Table 
^ and allow in principle for an extraction of the whole combination of LECs. X is 
counted in units of 10~ 2 . The results collected in Table |2| exhibit very clearly the 



Table 2: Contributions of loops and electric counterterms to the branching ratio of 
Kl — > 7r + 7r~e + e~ for different cuts in q 2 and for the entire phase space, given as 
functions of the combination (Vf 4 (/i) - JVJ" 6 (/i) - 3(JVJ" 6 (/i) - N 17 )) =: X(»[l(r 2 ] of 
weak LECs. The error is due to the uncertainty of Lg. 



q 2 > (MeV 2 ) 


Loops+Counterterms BR [10 8 ] 


2 2 


0.87 + 0.46V + 0.06V 2 


10 2 


0.86 + 0.45V + 0.06V 2 


20 2 


0.84 + 0.44V + 0.06V 2 


30 2 


0.80 + 0.42V + 0.05V 2 


40 2 


0.75 + 0.39V + 0.05V 2 


60 2 


0.61 + 0.32V + 0.04V 2 


80 2 


0.46 + 0.23V + 0.03V 2 


100 2 


0.30 + 0.16V + 0.02V 2 


120 2 


0.18 + 0.09V + 0.01V 2 


180 2 


0.007 + 0.003V + 0.0004V 2 


entire p.s. 


0.87 ± 0.19 + (0.46 ± 0.03)V + 0.06V 2 
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entire contribution of the involved weak local counterterms to the electric order p 4 
branching ratios. 

With the help of Eq. (|23|), we extract for the combination {N[ A (m p ) — Nl 5 (m p )) 
a central value of —0.019; therefore, the next-to-leading-order electric form factors 
effectively depend only on (Nl e (fi) — Nu) =: x(p) and electric contributions of order 
p 4 to the branching ratios can also be expressed as functions of x, again counted in 
units of 10~ 2 . The derived branching ratios are listed in Table ^. The error in Table 
^ was estimated by taking into account the uncertainties of Lg and (-/V{" 4 — N[ 5 ). 

Table 3: Contributions of loops and electric counterterms to the branching ratio of 
Kl — > 7r + 7r~e + e~ for different cuts in q 2 and for the entire phase space, given as 
functions of the combination (iVf 6 — N±y) =: x [10 -2 ] of weak LECs. (A r { 4 (mp) — 
JVf 5 (m p )) = -0.019 was used. 



q 2 > (MeV 2 ) 


Loops+Counterterms BR [10 8 ] 


2 2 


0.22 - 0.68x + 0.54x 2 


10 2 


0.22 - 0.67x + 0.53x 2 


20 2 


0.21 - 0.66a; + 0.52x 2 


30 2 


0.20 - 0.63x + 0.49x 2 


40 2 


0.19 -0.59x + 0.46x 2 


60 2 


0.16 -0.48x + 0.37x 2 


80 2 


0.12-0.36x + 0.27x 2 


100 2 


0.08 - 0.24x + 0.18x 2 


120 2 


0.05-0.14x + 0.10x 2 


180 2 


0.002 - 0.005x + 0.004x 2 


entire p.s. 


0.22 ± 0.11 - (0.68 ± 0.16)x + 0.54x 2 



The numbers in Tables |^ and ^ cannot be compared immediately to the results 
in H , since the corresponding branching ratios were expressed as functions of a dif- 
ferent combination of LECs, wl |§- It turns out that wl is related to the used iVj 
through w L = 8tt 2 [-7Vi4 + N 15 + N 16 - N 17 }. 

It is obvious from Tables ^ and |3| that the electric 0(p 4 ) contributions are nearly 
insensitive to changes of the cut below ~ (40 MeV) 2 . This feature becomes also 
clear from inspection of Figure ||. 

Theory finally predicts as central value of BR(Kl — > tt + tt~ e + e~) over the entire 
phase space [21.2 + 12.8 + 0.87 + 0.46X + 0.06X 2 ] • 10~ 8 , with X = {N{ A - iVf 5 - 
3(iV{ 6 — iV"i7))[10 -2 ] . Comparison with the branching ratio obtained in Ref. Q, 
BR = [18 (magn.) + 13 (tree) + 0.4 (CR)] • 10~ 8 , shows that inclusion of the mag- 
netic form factor of Ref. []llj increases the magnetic BR considerably. Also the total 
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0(p 4 ) electric contribution of Table [2| changes the result to some extent. 

In the following, the obtained theoretical BR over the entire phase space from 
Table |3] will be compared to the most recent available data to extract values for 
(Nl 6 (m p ) — Nu). It should be mentioned that possible values of the related combi- 
nation wl of LECs were estimated in Ref. j|] by comparing with the then recent BR. 
However, a theoretical cut of q 2 = (2 MeV) 2 was applied and the energy dependence 
of the magnetic form factor (19) was not taken into account. 



Figure 6: Differential decay width y^~w ^ or ( N i6( m p)~ N ^) = -0.014; r := q 2 /m 2 . 
The dotted line is the tree-level contribution, the dashed line refers to loops and 
counterterms, the dot-dashed line refers to the magnetic part. The thick long-dashed 
line is the sum. Tk l is the total width of the K^. For cuts with q 2 > (130 MeV) 2 
the differential decay width is dominated by the contributions of the electric 0(p 4 ) 
amplitude, i.e. loops and electric counterterms. 



I focus on the data of the KTeV collaboration, but for completeness one should 
mention that a Japanese group obtained a BR of [4.4 ± 1.3(stat.) ± 0.5(syst.)] • 10~ 7 
H, based on 13 events, and that the NA48 experiment at CERN recently reported 
a preliminary BR of (3.1 ± 0.3) • 10~ 7 [0, gj. In the last years, the KTeV result for 
the branching ratio was subject to numerous analyses and the errors improved quite 
a lot. The first published BR was based on a sample of 46 events and it was found 
to be [3.2 ± 0.6(stat.) ± 0.4(syst.)] • 10~ 7 [§. A new analysis based on the full 1997 
data set reported a BR of [3.32 ± 0.14(stat.) ± 0.28(syst.)j • 10~ 7 with a much better 



statistical error [12|. I am going to use the latest available (preliminary) numbers 
which were again obtained from the 1997 data set by considering the parametriza- 
tion in © fjj], IH : BR = t 3 - 63 ± 0.11(stat.) ± 0.14(syst.)] • 10" 7 . 
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It is clear that it is not possible to determine unambiguously the value of x only 
by comparison with the experiment. The two possible values of x are 



{N{ & (m p ) - N 17 )x = Xl = (2.7 ± 3.6) • 10~ 2 , 

(Nl 6 (m p ) - N 17 ) 2 = x 2 = (-1.4 ±3.6) -10-2. (38) 

The large error is mostly (~ 80%) due to the uncertainty of the magnetic BR. There 
is even a small overlap of the two ranges of x% and x 2 - Moreover, it should be 
stressed that comparison with the BR over the entire phase space is not the best 
possibility to extract values for the LECs, since the BR over the whole phase space 
is dominated by the tree level and the magnetic amplitude. In addition, Figure 
^ suggests that one could extract a value to a better precision for much higher 
cuts in q 2 , but this is not possible at the moment. On the other hand, such an 
extraction would suffer from smaller statistics. Nevertheless, the central values in 



(38) are very different and one can appeal to models of weak counterterm couplings 
to distinguish between the two solutions. The Weak Deformation Model (WDM) and 
the Factorization Model (FM) [2^] make predictions about the involved LECs, but 



apart from a free parameter of the FM, iY{* 4 and N[ e depend in both models on the 



contact term coupling H\ from the strong counterterm Lagrangian [19]. Therefore, 
it is necessary to compare (N{± — N{ 5 ) — (iV{ 6 — N% 7 ) to the experimental values, since 
in this combination H\ drops out. The comparison of prediction and experiment 



Table 4: Comparison of model predictions (WDM, FM) and the two possible values 
extracted from data for the combination (iVf 4 — N[ 5 ) — (JVfg — Ni 7 ) at m p . 



Model 


Pred. 


(NU-N[ 5 )- Xl 


(NU-N{ 5 )-x 2 


WDM 


-0.004 






FM 


-0.007&/ 


-0.046 ± 0.036 


-0.005 ± 0.036 



is given in Table |j. kf parametrizes the factorization hypothesis and is expected 
to be of O(l). Comparison with the results in Table [| gives kf\ ~ 6.4 ± 5.0 and 
kf2 — 0.7 ±5.0, respectively. First of all, it is remarkable to find the central value of 
X2 to be in such good agreement with the predictions of the two models. Secondly, 
the errors are big enough to dampen too much enthusiasm, but in any case the 
solution (iV{ 6 (?7ip) — N\ 7 ) = —0.014 is clearly favoured by both models. Moreover, 
both models and also the results of [33] suggest that Nn vanishes individually, so 



one can even go one step further and assume that 

N{ 6 (m p ) = (-1.4 ± 3.6) • 10~ 2 , iV"i 7 = 0. (39) 
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Additional support comes from the null measurement of interference between electric 
direct emission and bremsstrahlung in the K + — > 7r + 7r°7 amplitude in Ref. [17], since 



this indicates that the combination (-/V{* 4 — iV{ 5 — N^g — Nn) ~ or very small |23|, E5|. 
Assumption (j39|), however, is only true for a certain class of models and one should 
take into account other approaches, too, e.g. the modified FM (FMV) approach 



[34, 35] which was originally introduced to estimate 0(p ) corrections to radiative 
kaon decays. This model can also be used to parametrize LECs and it predicts in 
general a Nn different from zero. In any case, solution X2 is also supported by the 



results in 35]. We will consider the FMV more closely in the next section. 

If one's trust in the models used above were big enough, one could even use 
(|39|) to calculate the contact term coupling Hi of Ref. |n| and, as a consequence, 
calculate iV{ 4 and N[ 5 , but this does not seem to make much sense. In any case, 
assumption (|39|) serves as a solid starting point for the analysis of the K + decay. 



4.3 Numerical analysis of K + — > ir + n°e + e 

As in the previous case, there is no interference between magnetic and electric parts 
of the amplitude. Because of the absence of a symmetry relation between the electric 
form factors as in the Kl decay, however, this time there is interference between the 
tree- level amplitude and loops and electric counterterms. 

For the purely magnetic part of the branching ratio, the value of (|o|) is used. 
The results are collected in the second column of Table [|. As in Sect. 4.2, I only 
quote the error associated with \A±\ for the branching ratio over the entire phase 
space. 

In the following, I present the individual BRs due to electric form factors of 
lowest and next-to-leading order as well as the total electric branching ratios, which 
allows for an extraction of the interference contribution. The BR due to the lowest 



order is generated by the tree-level form factors given in (28). They produce a 
branching ratio that is much larger than that of Kl — > ir + Tr~e + e~; it is given in 
the third column of Table || As already mentioned, the tree-level value of G27 is 
used for the numerical analysis and this clearly introduces an intrinsic uncertainty 
in the predictions. The error associated with the choice of G$ is very small for the 
tree level. For loops and counterterm contributions we choose again the canonical 
\Gg\. For completeness, I quote the error arising from numerics because of the l/q 4 
behaviour of the tree-level amplitude. 

The 0(p 4 ) form factors contain the combinations of LECs given in (32). Here, I 



use the assumption of (|39| ) and thus express the derived branching ratios as functions 
of (JVJiO) + 2iVf 5 (//)) =: z(n) (in units of 1(T 2 ). The results are given in Table 
^. The total electric contributions form the tree level, loops and counterterms are 
collected in Table |?|. From this analysis it is clear that it will be very difficult to 
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Table 5: Magnetic and tree- level contributions to the branching ratio of K + — > 
7r + 7r e + e _ for different cuts in q 2 and for the entire phase space. The error of the 
magnetic part is due to experimental uncertainties of \A^\, whereas the error of the 
tree-level result comes from numerics; additionally, there is an intrinsic uncertainty 
because of the G27 coupling. 



q 2 > (MeV 2 ) 


Magnetic BR [1(T 8 ] 


Tree-level BR [10" 8 ] 


2 2 


5.33 


254.20 


10 2 


2.84 


74.33 


20 2 


1.80 


32.51 


30 2 


1.23 


17.35 


40 2 


0.86 


10.04 


60 2 


0.42 


3.75 


80 2 


0.19 


1.46 


100 2 


0.083 


0.56 


120 2 


0.031 


0.20 


180 2 


0.0002 


0.002 


entire p.s. 


6.14 ±1.3 


330 ± 15 



isolate the electric 0(p 4 ) corrections to branching ratios with small or no cuts in q 2 . 
Again, the importance of the last step from a cut of 4 MeV 2 to no cut at all should 
be mentioned. Comparison with the previous Kl decay shows that the tree-level 
contribution, although it is suppressed by isospin symmetry, dominates the BR and 
that it is much more important than for the Kl decay, where the tree level was 
e-suppressed. 

Finally, using only the central values of input quantities and applying assumption 
(|39|), Nl e (m p ) = —0.014 and N\j = 0, the central value of the branching ratio for 
K+ — > 7T + 7r°e + e over the entire phase space is predicted to be [6 + 378 + 0.27z + 
0.004z 2 ] • 10~ 8 . 



4.4 Dependence on counterterm models 

In Sect. 4.2 it was claimed that also the modified Factorization Model (FMV) of 



Refs. @, H favours (N[ 6 (m f 



N- 
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-0.014. We will now clarify why this is 



so. In Refs. []34, 35], the authors introduced a different approach compared to the 
one used in [22] to incorporate interactions between pseudoscalars and vector and 
axial-vector resonances. Factorization, however, was still an important ingredient. 
Originally used to estimate order p e corrections to kaon decays, their framework 
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Table 6: Electric 0(p 4 ) contribution to the branching ratio of K + — ► 7r + 7r°e + e for 
cuts in q 2 and for the entire phase space considering the error of (|39j). (N[ 4 +2N[ 5 ) := 
z [HT 2 ]. 



q 2 > (MeV 2 ) 


Loops+Counterterms BR [10 9 ] 


2 2 


46 + 0.51z + 0.042z 2 


10 2 


2.52 + 0.48z + 0.042z 2 


20 2 


1.68 + 0.46z + 0.041z 2 


30 2 


1.43 + 0.43,2 + 0.039z 2 


40 2 


1.25 + 0.40z + 0.036z 2 


60 2 


0.94 + 0.32z + 0.0302 2 


80 2 


0.67 + 0.24z + 0.022z 2 


100 2 


0.43 + 0.16,2 + 0.015z 2 


120 2 


0.25 + 0.09z + 0.0092 2 


180 2 


0.01 + 0.004z + 0.00042 2 


entire p.s. 


122 ± 134 + (0.56 ± 0.27)z + 0.043z 2 



was also extended to parametrize combinations of weak LECs in terms of two pos- 



itive 0(1) parameters, r/y and r/A- According to [35|, however, one finds that the 
weak LECs we are interested in, N{ 4 , N[ 5 , N[ e and Nu, do not depend on the 
factorization hypothesis. Hence, one should consider the parametrizations of combi- 
nations of these LECs in terms of rjy and t\a as model independent. Of course, these 
parametrizations of combinations of LECs still depend on the formalism applied to 
incorporate vector and axial-vector resonances |35| Q 

In the FMV, the relations (iVf 4 - N{ 5 ) = -0.020^ + 0.004r/ A and (iVf 4 - iVf 5 ) - 
3(N'[q — Niy) = — 0.004?7y + 0.018?^ hold. Comparison with the two possible val- 
ues of (-/VJg — Nn) in (^) gives the following results for the two parameters: for 
{N[ 6 {m p ) - N 17 ) = 0.027, one finds 

77^= -0.2 ±1.1 and r) Al = -5.6 ± 5.8. (40) 

Using (A^{" 6 (mp) — ^17) = —0.014, one calculates for the FMV parameters 

7/y 2 = 1.3±1.1 and t]a 2 = 1.6 ± 5.8. (41) 

Of course, the errors are very big, but even then the second pair of values in (|4l|) fits 



much better than the values of (40). Using the parametrization (iV{" 4 — iVJ 5 ) — 3(iV{ 6 + 



thank J. Portoles for useful comments on this topic. 
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Table 7: Contributions of the total electric part of the amplitude, 0(p 2 ) and 0(p 4 ), 
to the branching ratio of K + — » 7r + 7r°e + e _ for cuts in q 2 and for the entire phase 
space with consideration of the error of (|39|). (N[ 4 + 2A r { 5 ) := z [10~ 2 ]. 



q 2 > (MeV 2 ) 


Electric 0{p 2 ) + 0(p 4 ) BR [1(T 9 ] 


2 2 


2745 + 2.69z + 0.042z 2 


10 2 


785.6 + 2.57z + 0.042z 2 


20 2 


354.3 + 2.29z + 0.041z 2 


30 2 


193.8 + 2.09z + 0.039z 2 


40 2 


115.6 + 1.83z + 0.036z 2 


60 2 


46.3 + 1.32z + 0.030z 2 


80 2 


19.8 + 0.88z + 0.022z 2 


100 2 


8.40 + 0.52z + 0.015z 2 


120 2 


3.34 + 0.27z + 0.009z 2 


180 2 


0.06 + 0.008z + 0.0004z 2 


entire p.s. 


3783 ± 350 + (2.74 ± 0.82)z + 0.043z 2 



N 17 ) = 0.05??y - OMrjA, one determines N 17 = -0.009r/y + 0.0097^ = 0.4 • 10" 



Thus, we rather find for N[ e and Nn with the values of (41) 



ATeK) = (-1.0±4.6)-10" 2 , 

JV17 = (0.4 ± 4.6) • 10" 2 (42) 

In fact, this result is not too different from the hypothesis in (|3S|). 

Despite the big uncertainties of the values of the parameters rjy and t\a-, one 
nevertheless can use the central values to calculate other LECs, especially N[± and 
N[ 5 . Probably this is equally daring as the option of calculating H\ , but if one truly 



'believed' in the FMV with the results of (41), one could make a parameter- free 
prediction for the electric part of the branching ratio. According to [35], we have 
2iVf 4 + iVf 5 = 0.087/y. Using the experimental result (iVJ 4 - N{ 5 ) = -0.019, one 
calculates: 

N[ A {m p ) = 2.8 -HT 2 , 

N{ b {m p ) = 4.7 -10" 2 . (43) 

The results for the electric branching ratios obtained with these values for the coun- 
terterm couplings are listed in Table ||. According to Tables ^ and |?], one finds that 
0(p A ) corrections become more important for higher cuts in q 2 and that branching 
ratios for lower cuts are dominated by the tree level. E.g. for a cut of (10 MeV) 2 , 
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Table 8: Total electric branching ratio of K + — ► 7r + 7r°e + e _ for cuts in q 2 and for 
the entire phase space relying on the values of the LECs calculated in (|42]) and ( [4S|) . 
No errors are taken into account. 



g 2 > (MeV 2 ) Electric C(p 2 ) + 0(p 4 ) BR [10~ 9 ] 



2 2 


2923 


10 2 


873 


20 2 


431 


30 2 


253 


40 2 


164 


60 2 


78 


80 2 


39 


100 2 


20 


120 2 


9 


180 2 


0.24 


entire p.s. 


3683 



the tree-level BR is modified by 0(p 4 ) corrections and by the interference between 
the two electric contributions by about 17%, whereas for a cut of (80 MeV) 2 the 
result is increased by roughly 160%. It is clear that the interference also gives rise 
to an important part of the enhancement of the BR. Finally, one should note that 
the estimated couplings in (|42|) and (|43|) are in a range where one could expect them 
but it is also true that the uncertainties involved are too big to make a more precise 
statement about the couplings and the K + decay width. 

5 Conclusions 

I considered the non-leptonic decays Kl — > 7r + 7r _ e + e _ and K + — > 7r + 7r°e + e _ within 
the framework of Chiral Perturbation Theory. First of all, the amplitudes of the de- 
cays have been given up to order p 4 in a very explicit way and a consistency check on 
parts of the weak counterterm Lagrangian of CHPT was performed: all divergences 
are properly removed. 

The main reason to focus on Kl — > ir + Tr~e + e~ in this paper is provided 
by the possibility of extracting the combination (iVJ 6 (/z) — N\j) of weak LECs 
from experimental results. The latest value of the preliminary branching ratio, 
BR(K L -> Tr+vr-e+e") = (3.63 ± 0.11 ± 0.14) • 10~ 7 g§ and the values of the 
parameters of the magnetic form factor [pH[ , both obtained by the KTeV collabora- 
tion, were used for the numerical analysis of the decay. 
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The introduction of an energy dependent magnetic form factor yields an impor- 
tant correction to the older calculations in [|], ||, |3| , since it increases the magnetic 
contribution to the branching ratio considerably. The (preliminary) value for the 
branching ratio and the parameters of the magnetic form factor were obtained from 



the analysis of the data set of 1997 which contains more than 1800 events [11 
Comparison with earlier experimental results shows that the errors of the measured 
quantities became quite smaller due to the better statistics but the uncertainties are 
still too big to make precise predictions. 

Comparison with experiment yields two possible values of the LEC combination 
(N[ G (m p ) — Nn), thus one has to consult theoretical approaches about weak coun- 
terterm couplings to distinguish between the possible solutions. All models that 
have been considered (Weak Deformation Model WDM [22], Factorization Model 
FM Itffi l, modified Factorization Model FMV [34, |3^[) prefer the same value of 
(—1.4 ± 3.6) • 10~ 2 . Of course, the error, which is mainly due to the experimental 
uncertainties of the two parameters of the magnetic amplitude, is quite large, but 
nevertheless the obtained result is reasonable compared to (N[ 4 (m p ) — N[ 5 (m p )) = 
— 1.9 • 10 -2 and it rests upon a firm theoretical ground. 

On the other hand, the central value of (iVfg — Nn) is almost in perfect agree- 
ment with the FM and WDM predictions. One also derives central values for the 
two parameters of the FMV that are in good agreement with the expectations. 

Since 1997, much more data have been collected by the KTeV group and there- 
fore one can hope that a new analysis of the much bigger set of events can reduce the 
experimental uncertainties. As already pointed out in Sect. 4.2, it should also be 
possible to extract the value of (-/V{ 6 — Nn) to a better precision by comparing the 
theoretical results with branching ratios for higher cuts in q 2 (e.g. ~ (40 MeV) 2 ), 
since the contributions from loops and counterterms become much more important 
for higher cuts than for the entire phase space. 

To be able to make a useful prediction for the K + decay, one has to rely on addi- 
tional theoretical assumptions. First, I followed the predictions of the FM and WDM 
and assumed that the extracted value —1.4- 10~ 2 is produced solely by N[ 6 and that 
Nyj = 0. It is therefore possible to express the branching ratio for K + — ► 7r + 7r°e + e _ 
as a function of (N[ A + 2iV{ 5 ). 

Contrary to the decay of the Kl, there exists an interference between the tree- 
level amplitude and the electric 0(p 4 ) amplitude. In general, it is found that the 
branching ratio due to the electric part of the decay amplitude clearly dominates 
over the magnetic contributions. Moreover, it is the tree level that produces by far 
the most important contributions to the BR over the entire phase space as well as 
for a wide range of cuts in q 2 . 

An extraction of the combination (Nl 4 + 2N[ 5 ) from the branching ratio over the 
entire phase space is almost impossible, but according to the discussion in Sect. 4.4, 
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with a scan of the q 2 spectrum it is more likely to extract values for (iVf 4 + 2iVf 5 ), 
especially for cuts larger than ~ (60 MeV) 2 . Obviously the experimental error of 
the magnetic part of the amplitude and the error of the combination (-/VJ 6 — Nn) 
will not make it easier to extract a reasonable value, but hopefully new results from 
KTeV (and from CERN) also help to improve the predictive power of this analysis. 

Whereas the analysis summarized so far was based on conservative assumptions, 
I also speculated about extracting values for iVj^, N'( 5 , N'[ e and -/V17. Referring to 
the FMV, the two parameters of the model were estimated using the available data 
and the extracted value of (A^{ 6 — iV 17 ). The central values of these parameters were 
further used to estimate the central values of the four low-energy couplings N{ 4 , 
N[ 5 , iVJg and Nn and to make a 'prediction' of the K + branching ratio without any 
free parameter. Although the errors are big, the obtained values for the LECs seem 
to be reasonable. 
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Appendix 

A Loop functions 

All loop integrals in this work can be reduced to a basis of three scalar integrals: 



where I introduced the abbreviations D\ = [k 2 — m 2 ], D2 = [k 2 — m 2 ][(k — q) 2 — M 2 ] 
and D 3 = [k 2 — m 2 ][(k — q) 2 — M 2 ][(k— p) 2 — M 2 }. Indexed loop functions are defined 




(44) 



27 



through the following relations: 
d d k k„ 



(2ir) d D 2 
d Jc 

d }\i kjjky 
(2vr) d D 2 



^£i(g 2 ,m 2 ,M 2 ), 
iq^Ci (q 2 ,p 2 ,qp,m 2 ,M 2 ) + ip^C 2 (q 2 ,p 2 , qp,m 2 , M 2 ) 
ig^B 00 (q 2 , m 2 ,M 2 ) + iq^Bniq 2 ,m 2 , M 2 ), 



(45) 



(2vr) 2 D 3 ~ i9 ^ Coo ^ q > p >QP> m ' M )+ i QnQvCn(q ,p ,qp,m ,M ) 

+i (Qij.Pu + QuPfj, )Ci 2 (q 2 ,p 2 ,qp,m 2 ,M 2 ) + ip^p u C 22 (q 2 ,p 2 ,qp,m 2 ,M 2 ). 

They can be given explicitly in terms of (^4|). Divergences arise through the scalar 
loop functions A and B in (fl4|). The divergent parts of these functions are isolated 
in expressions similar to (||). 



A(m 2 )\ div = -2m 2 A(ji), B(q 2 ,m 2 ,M 2 ) 
/~ 4 r 1 



div 



-2A0), 



16vr 2 



1 



d-4 



(ln(4vr) + 1 - j E ) 



(46) 



Apart from Coo, an C-like functions are finite. 

B K i — > 7r + 7r _ 7* form factor 

The electric form factor gets contributions from topologies 2 and 3 collected in 
J-'iq', the first part of is due to the loop particles (K®, ir^), the second part arises 
from (rjs, K^), and the very last line comes from (ir°, K^). 



-r 16 



2A(m|) - ^(q 2 + 2q Pl )B{q 2 ,ml,ml) + 



F 

2ml(q 2 + 2q Pl - m\) + m\{m\ - q 2 - 2qp,)\ (B (m» , m 2 K , m 

1 



-£(GPi + 9) 2 ,™?oO) + 



q 2 + 2qpi 



m 2 K (ml + pip 2 + QP2 - q 2 



-2qp x ) - 2ml(pip 2 + qp 2 ) Bi(ml,m 2 K ,rnl) + 



1 



2(g 2 + 2qp 1 ) 



m 2 K (-2ml 



-q 2 - 2qp\ - 2p\p 2 - 2qp 2 ) + m 2 r (2q 2 + Aqpi + Apip 2 + 4qp 2 ) + 2q 4 
+q 2 (8qpi + 2pip 2 + 2qp 2 ) + 8^pf + kp\p 2 qp\ + 4qp 2 qp 1 

B 1 (( Pl +q) 2 ,m 2 K ,m 2 7r )+4B 00 (q 2 ,ml,m 2 n ) + 



2 ...2 . ,2, , ,p ,.,:> _2 , 2(pip 2 ^ 1 qp 2 ) 

2 + 2qp\ 
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B 00 (ml,m 2 K ,ml) - B 00 {{pi + Q) 2 ,m 2 K ,ml) +m^n(m^m|,mj) 

1 



-{m^ + q + 2qp 1 )B 11 ({p 1 + q) ,m K ,m n 

i 2 ) + 2ml (.V 2 + 2?pi)] C(ml, (pi + q) 2 

-4m 2 + m K ~ 2 9 2 - ^qpi) + m 2 (4pip 2 + 4q 2 + 8gpi 

1 



m Ki m K ~Q ~ 2 QPi 



2ml) + 2m" (g 2 + 2q Pl ) C(ml, ( Pl + q) 2 , m 2 K , m 2 ) + - [m 2 K (-4 Pl p 2 



Ci(m 2 , {pi + q) 2 ,m 2 K ,ml) + - m|(-4mj - 6gpi - 4pip 2 - 4gp 2 + m 2 K 



-2q 2 ) +4ml(q 2 + 2q Pl + Pl p 2 + qp 2 ) + 2qpi{q 2 + 2q Pl ) 
C 2 {m\, (pi + q) 2 ,m 2 K ,ml) + m 2 K C 00 {ml, {p 1 + g) 2 , m 2 K ,ml) - 2(m 2 
-m|-)Coo(m 2 r , (p 2 + g) 2 , m\,m\) + (2m 2 K pip 2 + m|-m 2 - 2m 2 pip 2 ) 
Cn(m 2 , (pi + g) 2 , m 2 K ,ml) + m\{qpi + 2gp 2 + 2m 2 + 4p x p 2 ) - 2m 2 (gp 2 

+2pip 2 ) Ci 2 (m 2 , (pi + q) 2 ,m 2 K , ml) + m 2 K {2pip 2 + 2gp 2 + m 2 + gpi) 
-2m 2 (pip 2 + qp 2 ) C 22 (m 2 , (pi + q) 2 ,m 2 K ,ml) 
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6 A <"*> - 5 



to 2 + 2m|- + 6(2gpi + <ro 2 + p x p 2 + q 2 ) B (q 2 , m 2 K , m 2 K ) 
m l{- 6 (PiP2 + qP2 + Q 2 + 2gpi) + 2m|- - 11m 2 ) + 3m 2 



l8(q 2 + 2q Pl ) 

{ml + 6q 2 + 12gpi + 2m|- + 6pip 2 + 6gp 2 ) (i?(m 2 , m 2 , m 2 K ) 



-B{{pi + qY,m^m 2 K )) + 

,2i r1 2 n „2 



9(g 2 + 2gp!) 

+m 2 (-llm 2 - 9q 2 - \8q Pl - 4m 2 K - 39 Pl p 2 - 39gp 2 ) + 6m 2 K {q 2 + 2q Pl ) 

1 



9m v {m n + pip 2 + qp 2 -q - 2qp 1 

2 



Bi{m'l, m 2 , m^) - 9( y + 2 gpi) i 3m ^ 2g|2 + 4,?Pl + 3piP2 + 3gP2 + 3m ^ 
+m 2 (-llm 2 - 4m| - 29g 2 - 58q Pl - 39 Pl p 2 - 39qp 2 ) - 4m\{q 2 + 2qp x ) 
— 72q Pl — 72q qpi — 18q — 30q p\p 2 — 30(7 QP2 - 60qpiPi P 2 - 60qp\qp 2 

B 1 {{p 1 +q) 2 ,m 2 ,m 2 K )+5B 00 {q 2 ,m 2 K ,m 2 K ) + — ^ _^ 2gp ^ ( m l + 3 PiP2 
+3qp 2 )B 00 {ml, m 2 , m\) - {ml + q 2 + 2gp x + 3p x p 2 + 3gp 2 ) 

Booiipi + 9) 2 , rn 2 K )] + 3 ( g 2 + 2gpi) m ^ m ^ + 3pm + ^ qP2 ^ 
B u {ml,m 2 , m 2 K ) - (m 2 (m 2 + 2g 2 + 4gpi + 3p x p 2 + 3<?p 2 ) + q 2 {q 2 + 4<?pi 



29 



m 



+3pip 2 + 3<?p 2 ) + 4qPi + 6qPi(qP2 + P1P2)) Budp! + q) ,m v ,m K ) 
1 

~18 

; 1 ■> 2 2 , 2 , ., „>2 ,,,2 ,..2 , 

18 

,2 rM ro™ « c„2 o™2 , _2n , 2 ro ,2 



l{ - 6(pip 2 + q 2 + 2gpi + gp 2 ) - 11m 2 + 2m|-) + 3m 2 (6pip 2 

+6qp 2 + 6<? 2 + 12q Pl + m 2 + 2m^)J C(m£, (pi + q) 2 , m 2 , m z K ) + i [3m; 
(-14m 2 - 24pip 2 - 18gpi - 6<?p 2 - 6<? 2 - 2m|- + m 2 ) + m 2 (24m 2 K 
+30qp 2 + G6q P i + 48pip 2 + 30<? 2 + 11m 2 ) + m 2 K (36q Pl + 12gp 2 + 12g 2 



1 



2m 2 K 



+48p x p 2 +4m^)JC , i(m 2 ,(pi + q) z ,m 2 ,m 2 K ) + y^|3m, 2 (m 2 

-14m 2 - 12g 2 - 28qpi - 24qp 2 - 2A Pl p 2 ) + m 2 (11m 2 + 24m 2 K + 36q 2 
+76qpi + 48c/p 2 + 48pip 2 ) + 4m 2 K (6q 2 + 14gpi + 12gp 2 + 12pip 2 + m K ) 
+24qpi(q 2 + qp 2 + pip 2 + 2#pi) C 2 (m 2 , (pi + q) 2 ,m 2 ,m 2 K ) + 2(m|- - m 2 ) 

C o(m 2 , (p 2 + q) 2 ,m 2 ,m 2 K ) + ^(m 2 + m 2 K - m 2 )C o(m 2 , (pi + q) 2 ,m 2 ,m 2 K ) 



+ 3 



m 2 (5m 2 + 3gpi + 9pip 2 ) - m 2 (m 2 + 2m|- + gpi + 3pip 2 ) - 2m|- 



(<7Pi + 3pip 2 ) Cn(m 2 , (pi + q) 2 ,m 2 ,m 2 K ) + - m 2 (10m 2 + 3^ 2 + llgpi 

+9gp 2 + 18pip 2 ) - m 2 (2m 2 + 4m| r + g 2 + 3gpi + 3c/p 2 + 6pip 2 ) - 2m 2 K (q 2 

+3qpi + 3<?p 2 + 6pip 2 ) Ci 2 (m 2 , (pi + q) 2 , m 2 ,m 2 K ) + i m 2 (5m 2 + 3g 2 

+8gpi + 9<?p 2 + 9pip 2 ) - m 2 (m 2 + 2m 2 K + q 2 + 2q Pl + 3gp 2 + 3p x p 2 ) 
-2m 2 K (q 2 + 2qp 1 + 3gp 2 + 3p x p 2 ) C 22 (m 2 , ( Pl + q) 2 ,m 2 , m 2 K ) 

+ 2{d-l) A{m2K) + 4^1) 5 ^ 2 ' m ^' m ^)^ 2 " 4 "4)}- 



(47) 



C K + -> 7r+7r°7* form factors JF+' 16 , JF+ Z 26 

The electric form factors and get contributions from topologies 2 and 3. 
Contributions with an rjs in the loop are collected in expressions ^15^5, contribu- 
tions with a pair of any kaon and pion are collected in J 7 ^^, respectively. The 
contributions with an r)8 read as: 



T +l _ 
•M5 — 



9(q 2 + 2q Pl )(q 2 + 2<?pi + 2gp 2 



m 2 (6<? 2 + 12q Pl + 22q>p 2 ) 



+m^(-3g - Qqpi - qp 2 ) + 6q (q + 4gp x + Aqp 2 + pip 2 ) + 12qp 1 (2qp 1 



30 



+4qp 2 +P1P2) + 24qp 2 (qp 2 + P1P2) A(m 2 K ) + — 

J 18(g 2 +2gpi)(^ + 2gpi + 2gp 2 ) 

3m 2 (-2m 2 + m|- - q 2 - 2q Pl - 2 Pl p 2 ) + m 2 (26m 2 - 15m 2 K + 13g 2 
+26gpi + I2qp 2 + 38pip 2 ) + m 2 K (m 2 K - q 2 - 2qp\ - 6qp 2 - 8pip 2 ) + 6P1P2 
(q 2 + 2qp ± + 2<ro 2 + 2pip 2 ) ^(ra 2 ,) + ^( llm ^ ~ 8m K + ®QP2 + 6pij> 2 ) 



B(q 2 ,m 2 K ,m 2 K ) + 



mlm 2 K (-12m 2 + 24m^- 



54(g 2 + 2g Pl ) (q 2 + 2q Pl + 2qp 2 ) 
-58q 2 - 116gpi - 60pip 2 ) + m 2 pip 2 (-12g 2 - 24gpi - 240m 2 - 144pip 2 
-144qp 2 ) + m 2 K pip 2 (30m 2 K - 60q 2 - 120q Pl - 72qp 2 - 12p\p 2 ) + m 4 
(-28g 2 - 56gpi - 144gp 2 ) + m 4 K (5q 2 + !0q Pl + 36gp 2 ) + m 2 (60g 4 + 240gp 2 
+240g 2 gpi + 120q 2 qp 2 + 240qpiqp 2 ) + m 2 K (-24q 4 - 96qp 2 - 96q 2 qpi 

-48q 2 qp 2 - 9§qpiqp 2 ) + 3m|- - 96m^ £(m 2 , m 2 ,m 2 K ) + lg( - g2 ^ 2gyi ) 
m 2 (— 3m 2 + 8m|- — 2m 2 — 6qp 2 — 6pip 2 ) + m 2 (33m 2 - 24m|- + 18<ro 2 

+I8P1P2)] B((pi + g) 2 , m 2 , - 18(g2 + ^ + 

-2gpi - 2gp 2 - 2pxp 2 ) + m 2 (22m 2 + 8m 2 K + 27g 2 + 54q Pl + 54gp 2 
+78p lP2 ) + 6m^(g 2 + 2g m + 2gp 2 )l B x {ra\ , m 2 , rr&) 



9m 2 (-2m 2 - g 2 



+ 



(g 2 + 2qpi)(q 2 + 2gpi + 2gp 2 



m 2 m| r (-2m^ - m^ + + 2gpi - 4qp 2 



-2p\p 2 ) + m 2 pip 2 (20m 2 + 6g 2 + 12gpi + 12gp 2 + \2p\p 2 ) + m 2 K p\p 2 

{-m 2 K + q 2 + 2qpx + 2gp 2 + 2pip 2 ) + 4m 4 (g 2 + 2qp\ + 3qp 2 + 2m 2 ) - m%qp 2 



Bi{ml,m 2 K ,m 2 ) 



1 



9(g 2 +2gpi) 



m 2 (-10m 2 + 4m|- - 5g 2 - lOgpi + 18gp 2 



+18pip 2 ) + 4m K {2q + 4gpi - 3gp 2 - 3pip 2 ) + 6g (g + 4qpi + gp 2 
+P1P2) + I2qpi(2qp 1 + qp 2 +pip 2 ) B 1 {{p 1 + q) 2 ,m 2 ,m 2 K ) 

+2B 00 (q 2 ,m 2 K , m\) + -^^±^^B 00 (m 2 , m 2 , m 2 ,) 
2 

~ (g 2 + 2gpi)(g 2 + 2gpi + 2gp 2 

+hpip 2 ) - 3m 2 K (qp 2 + pip 2 ) + 3pip 2 (q 2 + 2gp x + 2gp 2 + 2pip 2 ) 



2m 2 (2m 2 - m 2 ^ + g 2 + 2gpi + 3gp 2 



J B 00 (m 2 , m^, m 2 ) - - ^ [2m 2 + 6gp 2 + 6 Pl p 2 + 2g 2 + 4q Pl 
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#00 ( (pi + q ) , m,, , m K ) + 3(g2 + 2qpi + 2gp2) 5 uK. m ,.%) 
2m 2 



2m 2 (2m 2 - m 2 K + q 2 + 2q Pl + 3<?p 2 



(q 2 + 2qp 1 )(q 2 + 2<roi + 2gp 2 
+5pip 2 ) - 3m 2 K (qp 2 + pip 2 ) + 3pip 2 (g 2 + 2gpi + 2gp 2 + 2pi£> 2 ) 



Bn(m n , m K ,m v ) - - ^ [m n (m n + 2q + Aqp x + 3qp 2 + 3pip 2 ) 

+q 2 (q 2 + Aqpi + 3qp 2 + 3p\p 2 ) + qpi(^qpi + 6qP2 + 6P1P2) 

Bu((pi + q) 2 , m 2 , m 2 K ) + — (m 2 (-3m 2 - 2m 2 + 8m 2 K - 6qp 2 - 6pm) 

+m 2 (33m 2 - 2Am 2 K + 18qp 2 + 18 Pl p 2 )) C(m 2 n , ( Pl + q) 2 , m 2 , m\) 

— ^m 2 (— 3m 2 + 6m 2 + 2Am 2 K + 18qp\ — 18qp 2 + 36pip 2 ) + m 2 (A9m 2 

— 30m|- - 6qp\ + 30qp 2 + 12pip 2 ) + m 2 K (—lQm 2 K — \2qp\ + \2qp 2 

-Mftf2))Ci« (pi + g) 2 ,«) - ( m ^(- 3m ^ + 6m 2 + 24ml 

+18<? 2 + 48gpi + 36gp 2 + 36pip 2 ) + m 2 (49m 2 - 30m|- - 6<? 2 + 32qp x 
+ 12qp 2 + \2p\p 2 ) + m 2 K (-l6m 2 K - \2q 2 - 56qp\ - 2Aqp 2 - 24pip 2 ) 

+24qp 1 qp 2 + 24qp 1 p 1 p 2 ^C 2 (ml,(p l + q) 2 ,m 2 ,m 2 K ) - - (-5m 2 + m 2 

+2m 2 r )C , oo("i 2 r, (Pi + q) 2 ,rn 2 , m 2 K ) + (m 2 (-5m 2 - 3gpi - 9pip 2 ) + "i 2 

(m 2 + 2ml + qpi + 3pip 2 ) + m 2 K (2qp 1 + 6pip 2 )) 

C n (m 2 , (ft + q) 2 ,m 2 v ,m 2 K ) + (m 2 (-5m 2 - 3g 2 - 8q Pl - 9qp 2 - 9 Pl p 2 ) 
+m 2 (m 2 + 2m|- + q 2 + 2gpi + 3qp 2 + 3pip 2 ) + m 2 K (2q 2 + 4gpi + 6gp 2 
+6 Pl p 2 ))c 22 (m 2 , ^ + g) 2 ,m 2 ,m|) + (m 2 (-10m 2 - Sg 2 - 11^ 
-9gp 2 - I8P1P2) + m^(2mj + 4m| + c/ 2 + 3<toi + 3gp 2 + 6pip 2 ) 
+m 2 K (2q 2 + 6gpi + 6qp 2 + 12pip 2 )) Ci 2 (m 2 , (pi + q) 2 ,m 2 ,m 2 K ) 

—ieGs 



I8(q 2 + 2gpi + 2gp 2 ) 



-3c/ 2 - 6croi - 6gp 2 + 12pip 2 + 5m| 



+ 10m 2 )^(m 2 f ) + —^————^54 Pl p 2 + 62m 2 - 29m 2 K )A(m 2 ) 
1 



54(q 2 + 2gpi + 2qp 2 ) 



m 2 (— 30g 2 — 60gpi — 60gp 2 — 60pip 2 — 44m 2 
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+58m 2 K ) + m 2 K (\2q 2 + 24q Pl + 24qp 2 + 24 Pl p 2 + 13m^)J B(m 2 n , m 2 , m 2 K ) 
f m 2 (6q 2 + 12qpi + 12gp 2 + 39pip 2 + HmJ + 4m|- 



9(q 2 + 2g Pl + 2gp 2 ) 



1 



-9m,) - 9piP2mjBi(m ff ,m„,m x ) + 2(g2 - — - ^ [m 7r (2m 7r 

+6pip2 + 9m|- - 6m 2 ) + m 2 K (-2m 2 K + 6pip 2 + 3m 2 ) - 6pip 2 m 2 
^(m^m^m 2 ) + 3B 00 (q 2 ,m 2 K ,m 2 K ) - 2B m ((pi + q) 2 ,m 2 ,m 2 K ) 



+ 



4(m 2 + 3pip 2 ) 



3(g 2 + 2gpi+2gff 2 ) jgoo(m7r ' m "' + 2jBoo((pi + ^ ' "V m ^ 
2 (2m 2 + 3pip 2 ) 2 2 2 Am 2 {m 2 + 3 Pl p 2 ) 

ql + 2 qPl + 2qp 2 B00{m ^ mK > m r>> + 3(9* + 2q Pl + 2qp 2 ) 



-nil 2 2 \ 

B n (m 7T ,m v ,m K ) - 



2m 2 (2m 2 + 3pip 2 ) 



q 2 + 2qpi + 2qp 2 
Coo(ml, (pi + g) 2 ,m 2 ,m^) j. 



B n (ml, m 2 K , m 2 ) + 2(m\ - ml) 

(49) 



The remaining contributions from (K, tt) pairs are collected in the form factors 
J 7 ^ and T^. The first parts refer to (K°, vr°) contributions, then (vr°, K~), (K°, tt - ) 
and (K + ,ir + ) contributions are presented. 



^16 



—ieGs f 1 r 1 / 2 / 2 2 

— ^ — 4 i—^ n—^ r \2mJ2m„+q +2qp\ 

F \(q 2 + 2qp 1 )(q 2 + 2qp 1 +2qp 2 )l 2 V - V - 

+2<?p 2 + 4pip 2 ) + m 2 K (-m 2 K + q 2 + 2qp 1 - 2qp 2 ) + 2p 1 p 2 (q 2 

2 

+2q Pl + 2 9 p 2 + 2p 1 p 2 )).4(m 2 ) + ^(2 m ^(2m^ - 2m| 
+q 2 + 2gpi + 2pip 2 ) + m 2 K (m 2 K - q 2 - 2qp x - 2p 1 p 2 )^j 
B(ml,m 2 K , ml) - (mlm 2 K (-m 2 K + 2m 2 + q 2 + 2gpi + 4<?p 2 + 6pip 2 ) 
-2m 2 j>ij> 2 (2m 2 + q 2 + 2gpi + 2gp 2 + 2p x p 2 ) + m 2 K p x p 2 (-m 2 K + <? 2 
+2gpi + 2qp 2 + 2pip 2 ) - Am\qp 2 - m 4 K qp 2 ^B 1 (ml, m 2 K , ml) 

-2{p\p 2 {2ml - m 2 K + q 2 + 2<roi + 2c/p 2 + 2pip 2 ) + 2m^qp 2 - m 2 K qp 2 

I r> I 2 2 2 \ , 2 D / 2 2 2 
(-BOO ("V , ™/0 "V ) + "V 8 !! ("V , m Ki m ir 



3(^ + 2qpi)(q z + 2gpi + 2qp 2 ) L 
(-<? 2 - 2gpi + gp 2 ) + 2g 2 (g 2 + 4gp x + 4gp 2 + pip 2 ) + Ap x p 2 (qpi 
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+2qp 2 ) + 8{qp\ + qp\ + 2qp 1 qp 2 )\ [ A{m 2 K ) - — _ - 
j4(m ^ ) + 6(rf- i) (g 2 ~ 4m 2 K )B(q 2 ,m 2 K , m\) 

~ i 2 | o u o ro V f 277 ^^ 2 + 2 <?Pi + 3gp2) - m|-(g 2 

(g 2 + 2qp l )(q 2 + 2gpi + 2gp 2 ) L 

+2gpi) + 2<? 2 (g 2 + Aqpx + 4gp 2 + P1P2) + 4pip 2 (<?Pi + 2qp 2 ) + 8{qp\ 
+qp 2 + 2gp lW2 )] .4(m 2 ) + g2 + ^ + 2qm [ ~ («Vk - ml) 

-m% + m 2 m,|- - 2m 2 (^i + qp 2 ) + 2m 2 K (qp 1 + qp 2 ))B(ml,m 2 K ,ml 

1 



g 2 + 2qp 1 



+2 Pl p 2 ( m 2 K -rnl)B 1 (rnlrn 2 K ,ml 

-2ml(qp 2 + pip 2 ) + 2m 2 K (qp 2 + Pl p 2 f)B(( Pl + g) 2 , 
+2(m 2 K - ml)(qp 2 + p 1 p 2 )B 1 ((p 1 + q) 2 , m 2 K , ml) 



( - 2m£ + 2m>| 



2 2> 

m K ,m ni 



+2B 00 (q 2 ,m 2 7r ,ml) + 2m 2 (m 2 - m# + + P1P2) - 2m 2 K (qp 2 
+P1P2) C(ml, {pi + q) 2 ,m 2 K , m 2 ) - 2m 2 (to 2 - m 2 K + qp 2 ) 
-2m 2 K qp 2 Ci(m 2 , (p x + q) 2 , m 2 K , ml) - 2m 2 (to 2 - to|-) 

C 2 (to 2 , (px + q) 2 ,m 2 K ,ml) - 2 Pl p 2 (ml - m 2 K ) 
C n (ml, (pi + q) 2 ,m 2 K ,ml) - 2(to 2 - m 2 K )(qp 2 + pip 2 ) 
C 22 (ml, (pi + q) 2 ,m 2 K ,ml) - 2(to 2 - m 2 K )(qp 2 + 2p x p 2 ) 
Ci2(^, (pi + q) 2 ,m 2 K ,ml) 

2m 2 (4m 2 - 4to|> - g 2 - 2gpi 



3(<? 2 + 2qpi)(q 2 + 2gpi + 2c/p 2 ) 
+8c/p 2 + 12pip 2 ) + m 2 K (2m 2 K + q 2 + 2qp\ - 8qp 2 - \2p\p 2 ) + 2pip 2 
(g 2 + 2gpi + 8<?p 2 + 8pip 2 )) A(m 2 ) - m 2 r (2m 2 (2m 2 - 3to|- - 2g 2 
-4gpi + 2gp 2 + 4pip 2 ) + m 2 K (2m\ + q 2 + 2qp\ - 2qp 2 - Qp\p 2 ) 
+4pip 2 (-<? 2 - 2gpi + gp 2 +p 1 p 2 )^B(ml,m 2 K ,ml) + 2(4to£ 

+4to^ (-2to|- - g 2 - 2gpi + qp 2 + 2pip 2 ) + 3to^(to 2 + <?p 2 + pip 2 ) 
+2m 2 pip 2 (-7m 2 r - 2g 2 - 4gpi + 2gp 2 + 2pip 2 ) - 3m 2 K pip 2 {q 2 
+2gpi + 2gp 2 + 2pip 2 ) - 8mlm 2 K qp 2 ^B 1 (ml,m 2 K ,ml) + 4(ro 2 
(2m 2 - m\ + q 2 + 2qp x + 6gp 2 + 8pip 2 ) - 3m 2 K (qp 2 + pip 2 ) 
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+3pip 2 {q 2 + 2qp! + 2qp 2 + 2p 1 p 2 )j ^B 00 (ml, m 2 K , ml) + m 2 



Bu(m^,m K ,m w ) J J +2^00(9 ,m K ,m K )j. 
T2I is organized in the same way and reads as follows: 



(50) 



•^26 



-ieGa 



q 2 + 2qpx + 2qp 2 



(2m 2 n ~ 3rn 2 K + 2p l p 2 )A(ml) + 



m 



K 



(2m 2 - m 2 K )B(ml, m 2 K ,ml) - (2m\ + m% - 2mW K - m 2 K p W2 ) 
Bi(ml,m 2 K ,ml) - 2pip 2 ^B 00 (ml,m 2 K ,ml) +mj5n(mj,m^mj)] 



+ 



Q(q 2 + 2#pi + 2c/p 2 ) 



- Ami - 3m 2 K + q 2 + 2gpi + 2<?p 2 - 4pip 2 



+ 



2(d- 1) 
1 



A(m 2 K ) + 



4(d-l) 



{q 2 -Am 2 K )B(q 2 ,m 2 K ,m 2 K ) 



q 2 + 2gpi + 2qp 2 



(ml + m 2 K + 2p lP2 )A(ml) + m 2 K (m\ - ml) 



B(ml m 2 K ,ml) + 2 Pl p 2 (m 2 K - ^(mj, m 2 K , ml) 
+2Boo(q 2 ,ml,ml) + 2(m 2 K - m 2 )Coo(«4, (pi + q) 2 ,m 2 K ,ml) 

+ 3(q 2 + 2q 1 p 1 + 2q P2 ) ^ ~ ^ + + 

+Pip 2 )B(ml, m 2 K ,ml) - 2[ml(2ml - 3m 2 K + 2pip 2 ) + m 2 K (2m 2 K 

2 



-pm^iml^ml 



+ 



K ~ m if + QPi + P1P2) 



B((p 2 + q) 2 ,m 2 K , ml) + (m 2 K + q 2 + 2gpi + 2gp 2 + 2pip 2 ) 

4(3m 2 - 2m|- + pip 2 ) 



Bi((p 2 + q) ,m K ,mJ 
+mlB 11 (ml,m 2 K ,ml) 



3(q 2 + 2qpi + 2gp 2 



n / 2, 2 2\ 



1 r 



(2m 2 - 3m|- + 2qpi + 2pip 2 ) 



£(g 2 , ml, ml) + (3m|- + 2g 2 + 4q Pl + 4gp 2 + 4 Pl p 2 ) 
Bi((p 2 + q) 2 ,m 2 K ,ml)\ - -m 2 K (ml - m\ + q Pl + pip 2 ) 

C(ml, ( P2 + q) 2 ,m 2 K ,ml) + I [ 4 m 2 (m 2 - m| + q P i + P1P2) + < 

- + gpi + qp 2 + Pl p 2 ) 



+4qp 2 (-m 2 K + qp x + Pip 2 ) + m^J(7 2 (m 2 , (p 2 + q) 2 ,m 2 K ,ml) - -m 2 K 
C 00 (ml,(p 2 + q) 2 ,m 2 K ,m 2 r ) + m 2 Cii(m 2 , (p 2 + g) 2 , , m 2 ) 
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+ + q P 2)C 22 (ml, (p 2 + q)\ m* K , ml) + (2ml + QP2) 



Cuirnl, (p 2 + < 



(51) 
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